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■ Abstract 



We discuss the Palatini formulation of modified gravity including a Yukawa-like term. It is 
. shown that in this formulation, the Yukawa term offers an explanation for the current exponential 

accelerated expansion of the universe and reduces to the standard Friedmann cosmology in the 
appropriate limit. We then discuss the scalar-tensor formulation of the model as a metric theory 
and show that the Yukawa term predicts a power-law acceleration at late-times. The Newtonian 
^ , limit of the theory is also discussed in context of the Palatini formalism. 
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O^' 1 Introduction 
u 

The expansion of the universe is currently undergoing a period of acceleration which is directly 
measured from the light-curves of several hundred type la supernovae [1, 2] and independently from 
^ ', observations of the cosmic microwave background (CMB) by the WMAP satellite [3] and other CMB 

?H I experiments [4, 5]. However, the mechanism responsible for this acceleration is not well understood 

and many authors introduce a mysterious cosmic fluid, the so called dark energy, to explain this 
effect [6]. Recently, within the framework of higher order gravity theories [10], it has been shown that 
such an accelerated expansion could be the result of a modification to the Einstein-Hilbert action 
[7]. One such modification has been proposed in [8] where a term of the form R^^ was added to the 
usual Einstein-Hilbert action, whose origin can be traced back to the string/M-theory [9]. It was 
then shown that this term could give rise to accelerating solutions of the field equations without dark 
energy. 

Based on this modified action, we first use the Palatini variational formalism to derive our field 
equations. In the Palatini formulation, instead of varying the action with respect to the metric, one 
views the metric and connection as independent field variables and vary the action with respect to both 
independently. In the original Einstein-Hilbert action, this approach gives the same field equations as 
the metric variation. However when nonlinear terms are added to the action. Palatini formalism leads 
to different dynamical equations. The importance of this formalism lies in the fact that, contrary 
to the theory resulting from the metric variation which explains the accelerated expansion of the 
universe but is in conflict with solar system experiments, it is free of such conflicts. Interestingly, the 
Palatini formalism has been shown to be equivalent to a scalar-tensor type theory [14] in which the 
scalar field kinetic term is absent from the action [11]. These results are important and fundamental 
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for this formalism. Also in [16], Nojiri and Odintsov have shown that a combination of the and 
terms can drive both the current acceleration and inflation. The Palatini form of this theory is 
studied in [13] where the authors obtain the current acceleration of the cosmic which is compatible 
with observations. In [22] the authors propose a Ini? term to account for the current acceleration 
and also inflation. The Palatini formulation of Ini? gravity is studied in [15] where the acceleration of 
the cosmic in the limit of small curvature is predicted. For a more recent and comprehensive review, 
the reader should consult [17] for models with modified gravity and [18] for models with the Palatini 
approach. 

In this paper, we have used the Palatini formalism in a theory where a Yukawa type term has 
been added to the action. Since the term R~^, mentioned above, is similar to the Yukawa term 
presented below when R is small, we also expect this theory to have its origins in string/M-theory, at 
least for small curvatures. We have shown that the resulting theory could account for the accelerated 
expansion of the universe. We also study our model in the context of what is known as the scalar- 
tensor theory [14] where the metric formalism is used, leading to a power-law acceleration in the 
limit of small curvature. The Newtonian limit of the theory is investigated and shown to be correctly 
satisfied. 



2 Friedmann equation in Palatini formalism 

To begin with, a quick review of the Palatini formalism [12, 13] would be in order. Consider the 
action 

5 = 1 d^x^C{R) + I d^x^CM, (1) 
where k = 8nG and Cm is the Lagrangian density for matter. Variation with respect to g^i, gives 

where a prime denotes differentiation with respect to R and Tf^,y is the energy-momentum tensor in 
the matter frame. We assume that the universe is filled with dust and radiation, represented by 
energy densities pm and pr respectively, thus 

= {-Pm - Pr,Pr,Pr,Pr}, (3) 



with T = g^^Tfj^i, = —pm, on account of = Pr/3- As is well known, in the Palatini formulation the 

connection is not associated with Qf^^i, but with h^i, = C'g^p which is known from varying the action 
with respect to F^^. The Christoffel symbols associated with h^j^n are given by 



r^. = {lu]^ + ^ [^^\,du)C' - g,.g^"d^C'] , (4) 
where the subscript g indicates association with metric g^i,. The Ricci curvature tensor is given by 

R^^Aa) - lc'-\^c'v,c' + c'-'v^v,c' + l-C'-'g^^VaV^C, (5) 



R 



with the curvature scalar written as 

R = R{g) + SC'Wi.VC' - ^C'^V^.C^'C', (6) 

where R^uig) is the Ricci tensor associated with g^i, and R = g'^'^R^^- Contracting (2), one gets 

jC'R -1L = -kT. (7) 
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Equations (5) and (6) now define the Ricci tensor with respect to h^i,. To derive the Modified 
Priedmann (MF) equation we follow [13, 15] and consider the fiat Robertson- Walker metric for the 
evolution of the cosmos 

ds'^ = -dt'^ + a{tfd:x.'^. (8) 

This choice is considered to be consistent with observations [2]. From equations (8) and (5) we get 
the non-vanishing components of the Ricci tensor 

i?oo = 3^ - lc'-\doC'f + ^C'-'VoVoC (9) 

Rij = -(^ad + 2d^ + c'~^r°jdoC' + y/:'" VoVo/:'j Sij. (lo) 

Substituting equations (9) and (10) into field equations (2), we find 

where = | is the Hubble parameter with p and p being the total energy density and total pressure 
respectively. If we solve R in terms of T from equation (7) and substitute it into the expressions for 
C and OqC in the above equation, we will get the MF equation. 



3 Yukawa term 

To continue, let us consider an action of the form 

C{R) = -R - 4 exp(-ai?), (12) 
R 

where a,/? > and, as can clearly be seen, the Yukawa term dominates at small curvature. This 
feature is important since, as well shall see, it helps to explain the accelerated expansion together 
with the assumption that < 0. The contracted field equation (7) now reads 

f{R) ^ ^ + a exp(-«i?) + I exp(-ai?) = -kT/P = (13) 

where f{R) is a monotonically decreasing function and we have that limij_>_oo /(-R) — > +oo and 
lim^^Q- f{R) —>■ —GO. Thus R is uniquely determined for any value of npm/P = x through equation 
(13). Prom the conservation equation 

Pm + ^Hpm = 0, (14) 

we find 

3/3 {a^/R + 2a/ R^ + 2/^^) ( Kpm\ 
° - [a2 + 3(a/i? + l/i?2) + exp(ai?)//3] \ P ) ^ ' 

Let us now define = F{x,a, f3)H such that in the limit /3 — > 0, F{x,q:,P) = 0, since from (15) 
one can see that OqC vanishes in this limit. Noting that R = R{x) and substituting this into equation 
(11) we get the MF equation 

jj2 ^ i^Pm + 2Kpr - (3[R/(3- exp ( - gfi) /.R] ^^^^ 



-I I a/3exp(— ai?) , /3exp(— ai?) 
^ ^ R + W 



[6 + 3F(x, a, /3) (1 + l/2F(x, a, /?))] 



It can be seen from equations (13), (15) and (16) that when (3 — > or equivalently F[x,a,(3) = 0, 
the MF equation will reduce to the standard Friedmann equation 

^2 ^ l^jpm + Pr) ^^^^ 
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Let us now consider the evolution of the universe in vacuum, i.e. pm = and pr = 0. If one defines 
the parameter A; > 1 according to the following relation 

i?o = — Infc, (18) 
a 

where the subscript zero represents association with vacuum and substitutes this into the vacuum 
field equation f{R) = 0, one gets 

In A; 

a = ±^= (19) 

Vfc/3(3-lnfe) ^ ' 

Assuming a > 0, we find 

Rq = -^Jkf3{3 - Ink). (20) 
Substituting this into the the vacuum MF equation and setting a; = results in 

VP 



Hi = ^^k{3-lnk). (21) 

Note that equation (16) reduces to the above equation in vacuum. This equation tells us that when 
^ ~ ~ {lO'^^eV)"^ and k < e^, modified gravity with a Yukawa term points to an accelerated 
evolution of the universe. It is worth noting that the role of is similar to a cosmological constant 
in this theory. If the density of dust cannot be ignored, i.e. npm/P ^ 1, then this will be equivalent 
to /? being as small as the dust energy, so that F{x,a,P) ~ and from equation (13) we obtain 
R^ —Kpm- Therefore, in the Big Bang Nucleosynthesis (BBN) epoch, the MF equation reduces to 
the standard Friedmann equation. Consequently, the Yukawa term would no longer be effective there. 
However, for the present epoch the condition Kpm/P S> 1 breaks down and the Yukawa term starts 
to dominate. 



4 Scalar-tensor formulation 



In this section we study the scalar-tensor formulation of our model along the lines introduced in [11]. 
In this theory, the Palatini form of the action is shown to be equivalent to a scalar-tensor type theory 
from which the the scalar field kinetic energy is absent. This is achieved by introducing a conformal 
transformation in which the conformal factor is taken as an auxiliary scalar field. We show that the 
resulting field equations in the small curvature regime predict a power-law type acceleration for the 
universe. Let tis then first give a brief review of the formalism discussed in [14] . 
Introducing the auxiliary scalar field 0, one may write the action (1) as [14, 16] 



S[gij.„,(j),tp„ 



-g[C{4>) + {R- ^)C'{4>)] + Sm[g^,u, ^r, 



(22) 



where iprn represents the matter field. The equation of motion resulting from this action for the 
scalar field is ^ = i?, as long as C"{(j)) / 0, making action (22) classically equivalent to the action 
(1). This action is sometimes referred to as the Jordan (physical) frame action. Under the conformal 
transformation 

gijiu = e\/¥*3^^, (23) 

with 



action (22) can be written in the Einstein frame as [14, 16] 



'5[5,ii/,^,V'n 



/ 



d x\J —g 



+ 5„ 



exp (-^2^/3$ j g^^, ip, 



(24) 



(25) 
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Here, 5^,^ is the metric in the Einstein frame [14] and R is the scalar curvature associated with g^i,. 
It may be shown that the potential is given by 

This is the standard form of the scalar-tensor type theories mentioned above. It is worth noting that 
since we are working with metric connections here, the kinetic term in action (25) is present, contrary 
to the case when one works within the framework of the Palatini formalism. 

Lets us now concentrate on the action represented by (12) for which the potential can be written 

as 

VU) = ^ ^. (27) 

2.(l + a/?^ + ^f 

When 4> is small, equation (24) is approximated by ~ ^/P exp and in this limit potential 

(27) takes the following form 

F($) ^ 4 ^ -V^exp (-J^^) . (28) 

As will be shown below, this potential predicts a universe evolving with a power-law acceleration. 
Typical plots of V{(p) are given in figure 1. For a^/? = 0.1, the potential has a maximum at ~ 0.47. 
Therefore, if the universe starts from the maximum, it either asymptotically evolves to a de Sitter 
solution, or undergoes a power-law acceleration. The minimum of the potential occurs at ~ —0.67 
and thus if the universe starts at this point, it is in an anti-de Sitter phase where the curvature is 
negative. Therefore, when the absolute value of (f) becomes small, the universe evolves to a power-law 
acceleration at late-times and at (p — ^ —3.98, V{4>) is unbounded from above. For a^/? = 0.6 the 
behavior of the potential for ^ > is qualitatively the same as in the previous case. In the region 
(f) < 0, when cf) — > —1.99 the potential becomes unbounded from below and for cf) 2± —2.1 it has a 
maximum. 

To proceed further, let us consider evolution of the scale factor with time. From action (25) one 
gets 

= K{p^+Pm), (29) 

*" + 3;J*' + ^-il^p„ = 0, (30) 

where a prime denotes d/dt, and 

p$ = ^$'2 + y($). (31) 

These are cosmological equations of motion in the Einstein frame where w is the usual equation of 
state parameter and pm is the matter density with H being the Hubble parameter [8, 16]. We must 
now solve the system of equations (29) and (30). We first consider the case where pm = 0. When the 
potential is given by (28) and cf) = R is small, a solution is given by 

a{i) oc (32) 

and 

4 ~ 

$oc--lnt. (33) 
Here, t is the time coordinate in the Einstein frame, which is related to the coordinate t in the Jordan 

_1 / 2k j. ~ 

frame by e 2 V 3 dt = dt. As a result 

3p/3 = t, (34) 
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Figure 1: The scalar potential given by equation (27) for = 0.1, left and = 0.6, right. 



The power-law acceleration also occurs in the physical (Jordan) frame 



a{t) 



_ 1 
e 2 



(35) 



which is consistent with the result in [8]. Hence, at small curvatures, the cosmic acceleration is 
predicted by the Yukawa-like term. 

Although there is no need to introduce dark energy in modified gravity, a better understanding 
of the equation of state in this theory is afforded by following [21] and define the effective Equation 
of State (EOS). Written in this form, it has the advantage that any parameterizations done within 
the context of EOS for dark energy can be compared with that of modified gravity. Let us then start 
from equation (35) and find the Hubble parameter as a function of the redshift z as 



H{z) = 2Ho{l + z 



.1/2 



(36) 



where a^/a = 1 + z with oq and Hq being the values of the parameter at the present epoch. It is worth 
noting that equation (36) is the same as that derived from the standard Friedmann equation with 
w = —2/3. Now, we may write the Friedmann equation in a formal fashion which would encapsulate 
any modification to the standard Friedmann equation in the last term regardless of its nature [21], 
that is 



H^H^ = + zf + SH^Hl 

where = p/poo Poc = SHq/k. Also, defining the effective EOS, denoted by Weff{z), as 



^eff{z) 



1 dln5H^ 

-1 + 77- 



3dln{l + z)' 

we can calculate Weff{z) using equations (36), (37) and (38) with the result 

1 [4-3(l + z)2fi^] (l + z) 



Weffiz) = -1 + 



3 4(1 + z)- 0^(1 + z)3 ■ 



(37) 



(38) 



(39) 



Figure 2 shows variations of w^ff{z) for various values of 0^. Equation (39) shows that — 1 < 
Wf.ff{z) < —2/3, which is consistent with the observational data when the limits imposed by appro- 
priate constraints are taken into account. It is therefore relevant at this point to have a brief look at 
these constraints and sec how our results can be interpreted. Before doing so however, it should be 
mentioned that the model presented here is not expected to produce an effective phantom behavior, 
at least for small curvatures, since for small R it behaves like R"^ theories for which it is impossible 
to have phantom behavior, although an effective quintessence behavior is feasible [23]. 

The present observational data [24] suggest that our universe is dominated by a mysterious form 
of dark energy. As a result, the universe expansion is undergoing a period of acceleration. In terms of 
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Figure 2: Variation of Weff{z) for fl^ = 0.1 solid line, fim = 0.2 dotted line and fim = 0.3, dashed line. 



the constant EOS parametrization, observational data provided by SDSS indicate that this constant 
is close to —1. In other words, the accelerating universe could be either caused by the cosmological 
constant w = —1, or quintessence, —1 < w < —1/3 or phantom era, w < —1. These constraints 
have been verified by the Supernovae, WMAP and cluster abundance observations [25]. The above 
results, derived within the framework of a scalar-tensor type model, seems to favor the quintessence 
which predicts values for w in the range given above. However, in the Palatini formulation presented 
in section 3, we have obtained the cosmological constant represented by equation (21), which seems 
to be more in line with the cosmological constant type models. 



5 Newtonian limit 

The Newtonian limit of theories that predict cosmic acceleration, that is in the limit of weak field 
approximation, is naturally of interest and should be examined since such limits are expected to 
be compatible with the present cosmological data. It is known that the criteria for the correct 
Newtonian limit of C{R) theories is provided by the Dick's condition for fourth order theories [19], 
that is C" (R)\ii=jiQ = 0. The same condition has also been obtained in the Palatini formalism when 
the equations of motion are of the second order [20]. With this criteria in mind, inspection of our 
Lagrangian, equation (12), suggests that it does not satisfy Dick's condition. However, this should 
not be considered as a setback since with a small modification which does not affect our general 
results, we should be able to obtain the correct Newtonian limit. 
Let us then modify Lagrangian (12) and write it as 



-aR 



Now, substitution of the above Lagrangian in equation (7) leads to Rq = ^Ink. Thus for k = 
we have that C"{Ro) = 0. This means that for the choice of k made above, we recover the correct 
Newtonian limit. The interesting point is that if one applies the Palatini formalism to the above 
Lagrangian, the second term does not contribute in general and Lagrangian (40) behaves exactly the 
same as our original Lagrangian does so that the predictions for a universe with accelerated expansion 
are unaltered. 



6 Conclusions 

In this paper we have discussed the Palatini formulation of modified gravity with a Yukawa like term. 
This term may be used to explain the current exponential accelerated expansion and reduces to the 
standard Priedmann equation in vacuum. The scalar-tensor formulation of our model was discussed 
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and shown to predict a power-law acceleration at late times. We have shown that in both formulations 
the accelerated expansion of the cosmic may be accounted for. However, we expect the universe to 
pass through a matter dominated era before the accelerated expansion phase is reached [26]. 

The Newtonian limit of this theory was also examined. It was shown that in the Palatini for- 
mulation our modified gravity action with a Yukawa term is equivalent to Lagrangian (40) , with the 
latter satisfying Dick's condition for the correct Newtonian limit. 
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